1 BCO comparison of CASAL and Casal2 model configurations

This document compares the results of at least 2 CASAL model configurations (base and at least one sensitivity) and up to 6 Casal2 model configurations (3 BetaDiff and 3 ADOL-C).

The CASAL model sensitivity 1 has a smaller minimisation tolerance value than the CASAL base model (1e-7 vs. 2e-5).

The Casal2 ADOL-C and BetaDiff low tolerance models have a smaller tolerance value than the CASAL base model (1e-7 vs. 2e-5).

1.1 BCO model characteristics

The main characteristics of the Test Case BCO (ling) CASAL model are:

  • one stock, ages 3 - 20
  • one area
  • years 1900 - 2019, projection year 2020
  • one time step
  • two stock categories: male and female
  • natural mortality-at-age is 0.14 for males and 0.20 for females
  • sex-specific logistic maturity
  • sex-specific von Bertalanffy age-length relationship
  • one length-weight relationship (\(W = aL^b\))
  • Beverton-Holt stock-recruitment relationship, with steepness (h) 0.75 and \(\sigma_R\) 0.6
  • no ageing error
  • one survey with sex-specific age-based logistic selectivity
  • three fisheries, all with sex-specific length-based logistic selectivity and retained/discarded catches

Observation data include:

  • survey index
  • fishery CPUE index
  • fishery proportions-at-length

Parameters estimated include:

  • B0
  • male \(L_{\infty}\)
  • male age-length CV2
  • sex-specific survey logistic selectivity (a50 and ato95)
  • YCS (recruitment deviations)

1.2 R environment

## [1] "Mon Sep 27 19:59:50 2021"
## R version 3.6.3 (2020-02-29)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Debian GNU/Linux 11 (bullseye)
## 
## Matrix products: default
## BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.9.0
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.9.0
## 
## locale:
##  [1] LC_CTYPE=en_NZ.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_NZ.UTF-8        LC_COLLATE=en_NZ.UTF-8    
##  [5] LC_MONETARY=en_NZ.UTF-8    LC_MESSAGES=en_NZ.UTF-8   
##  [7] LC_PAPER=en_NZ.UTF-8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_NZ.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] Casal2_21.08   casal_2.30     devtools_2.4.2 usethis_2.0.1  rlist_0.4.6.2 
##  [6] ggthemes_4.2.4 gridExtra_2.3  coda_0.19-4    ggmcmc_1.5.1.1 ggplot2_3.3.5 
## [11] tidyr_1.1.3    huxtable_5.4.0 dplyr_1.0.7    plyr_1.8.6    
## 
## loaded via a namespace (and not attached):
##  [1] Rcpp_1.0.7         lattice_0.20-41    prettyunits_1.1.1  ps_1.6.0          
##  [5] assertthat_0.2.1   rprojroot_2.0.2    digest_0.6.27      utf8_1.2.2        
##  [9] R6_2.5.1           evaluate_0.14      pillar_1.6.2       rlang_0.4.11      
## [13] rstudioapi_0.13    data.table_1.14.0  callr_3.7.0        jquerylib_0.1.4   
## [17] rmarkdown_2.11     desc_1.3.0         stringr_1.4.0      munsell_0.5.0     
## [21] compiler_3.6.3     xfun_0.26          pkgconfig_2.0.3    pkgbuild_1.2.0    
## [25] htmltools_0.5.2    tidyselect_1.1.1   tibble_3.1.4       reshape_0.8.8     
## [29] fansi_0.5.0        crayon_1.4.1       withr_2.4.2        grid_3.6.3        
## [33] jsonlite_1.7.2     GGally_2.1.2       gtable_0.3.0       lifecycle_1.0.0   
## [37] DBI_1.1.1          magrittr_2.0.1     scales_1.1.1       cli_3.0.1         
## [41] stringi_1.7.4      cachem_1.0.6       remotes_2.4.0      fs_1.5.0          
## [45] testthat_3.0.4     bslib_0.3.0        ellipsis_0.3.2     generics_0.1.0    
## [49] vctrs_0.3.8        RColorBrewer_1.1-2 tools_3.6.3        glue_1.4.2        
## [53] purrr_0.3.4        processx_3.5.2     pkgload_1.2.2      fastmap_1.1.0     
## [57] yaml_2.2.1         colorspace_2.0-2   sessioninfo_1.1.1  memoise_2.0.0     
## [61] knitr_1.34         sass_0.4.0

1.3 CASAL and Casal2 model output

source('../../R-functions/report_read_in_CASAL_MPD_files.R')
source('../../R-functions/report_read_in_Casal2_MPD_files.R')

1.4 Tables

Tables of parameter estimates and objective function components for the CASAL and Casal2 model MPD results

CASAL parameter estimates
rownamesBase_ModelSensitivity_1Percent_Diff
initialization.B0       14170.0000014170.000000.000
recruitment.YCS1     1.000001.000000.000
recruitment.YCS2     1.000001.000000.000
recruitment.YCS3     1.000001.000000.000
recruitment.YCS4     1.000001.000000.000
recruitment.YCS5     1.000001.000000.000
recruitment.YCS6     1.000001.000000.000
recruitment.YCS7     1.000001.000000.000
recruitment.YCS8     1.000001.000000.000
recruitment.YCS9     1.000001.000000.000
recruitment.YCS10    1.000001.000000.000
recruitment.YCS11    1.000001.000000.000
recruitment.YCS12    1.000001.000000.000
recruitment.YCS13    1.000001.000000.000
recruitment.YCS14    1.000001.000000.000
recruitment.YCS15    1.000001.000000.000
recruitment.YCS16    1.000001.000000.000
recruitment.YCS17    1.000001.000000.000
recruitment.YCS18    1.000001.000000.000
recruitment.YCS19    1.000001.000000.000
recruitment.YCS20    1.000001.000000.000
recruitment.YCS21    1.000001.000000.000
recruitment.YCS22    1.000001.000000.000
recruitment.YCS23    1.000001.000000.000
recruitment.YCS24    1.000001.000000.000
recruitment.YCS25    1.000001.000000.000
recruitment.YCS26    1.000001.000000.000
recruitment.YCS27    1.000001.000000.000
recruitment.YCS28    1.000001.000000.000
recruitment.YCS29    1.000001.000000.000
recruitment.YCS30    1.000001.000000.000
recruitment.YCS31    1.000001.000000.000
recruitment.YCS32    1.000001.000000.000
recruitment.YCS33    1.000001.000000.000
recruitment.YCS34    1.000001.000000.000
recruitment.YCS35    1.000001.000000.000
recruitment.YCS36    1.000001.000000.000
recruitment.YCS37    1.000001.000000.000
recruitment.YCS38    1.000001.000000.000
recruitment.YCS39    1.000001.000000.000
recruitment.YCS40    1.000001.000000.000
recruitment.YCS41    1.000001.000000.000
recruitment.YCS42    1.000001.000000.000
recruitment.YCS43    1.000001.000000.000
recruitment.YCS44    1.000001.000000.000
recruitment.YCS45    1.000001.000000.000
recruitment.YCS46    1.000001.000000.000
recruitment.YCS47    1.000001.000000.000
recruitment.YCS48    1.000001.000000.000
recruitment.YCS49    1.000001.000000.000
recruitment.YCS50    1.000001.000000.000
recruitment.YCS51    1.000001.000000.000
recruitment.YCS52    1.000001.000000.000
recruitment.YCS53    1.000001.000000.000
recruitment.YCS54    1.000001.000000.000
recruitment.YCS55    1.000001.000000.000
recruitment.YCS56    1.000001.000000.000
recruitment.YCS57    1.000001.000000.000
recruitment.YCS58    1.000001.000000.000
recruitment.YCS59    1.000001.000000.000
recruitment.YCS60    1.000001.000000.000
recruitment.YCS61    1.000001.000000.000
recruitment.YCS62    1.000001.000000.000
recruitment.YCS63    1.000001.000000.000
recruitment.YCS64    1.000001.000000.000
recruitment.YCS65    1.000001.000000.000
recruitment.YCS66    1.000001.000000.000
recruitment.YCS67    1.000001.000000.000
recruitment.YCS68    1.000001.000000.000
recruitment.YCS69    1.000001.000000.000
recruitment.YCS70    1.000001.000000.000
recruitment.YCS71    1.000001.000000.000
recruitment.YCS72    1.000001.000000.000
recruitment.YCS73    1.000001.000000.000
recruitment.YCS74    1.000001.000000.000
recruitment.YCS75    1.000001.000000.000
recruitment.YCS76    1.000001.000000.000
recruitment.YCS77    1.000001.000000.000
recruitment.YCS78    1.000001.000000.000
recruitment.YCS79    1.000001.000000.000
recruitment.YCS80    1.000001.000000.000
recruitment.YCS81    1.000001.000000.000
recruitment.YCS82    1.000001.000000.000
recruitment.YCS83    1.000001.000000.000
recruitment.YCS84    0.234800.234800.000
recruitment.YCS85    0.231600.231600.000
recruitment.YCS86    0.232200.232200.000
recruitment.YCS87    0.241200.241200.000
recruitment.YCS88    0.265100.265100.000
recruitment.YCS89    0.317200.317200.000
recruitment.YCS90    0.422400.422400.000
recruitment.YCS91    0.613900.613900.000
recruitment.YCS92    0.923000.923000.000
recruitment.YCS93    1.429001.429000.000
recruitment.YCS94    1.480001.480000.000
recruitment.YCS95    1.009001.009000.000
recruitment.YCS96    0.846200.846200.000
recruitment.YCS97    0.855500.855500.000
recruitment.YCS98    1.087001.087000.000
recruitment.YCS99    1.530001.530000.000
recruitment.YCS100   1.027001.027000.000
recruitment.YCS101   0.799900.799900.000
recruitment.YCS102   0.820800.820800.000
recruitment.YCS103   0.631800.631800.000
recruitment.YCS104   0.805200.805200.000
recruitment.YCS105   0.745300.745300.000
recruitment.YCS106   0.463400.463400.000
recruitment.YCS107   1.075001.075000.000
recruitment.YCS108   0.899200.899200.000
recruitment.YCS109   0.385100.385100.000
recruitment.YCS110   0.912700.912700.000
recruitment.YCS111   0.709200.709200.000
recruitment.YCS112   0.829000.829000.000
recruitment.YCS113   0.423300.423300.000
recruitment.YCS114   1.278001.278000.000
recruitment.YCS115   0.688100.688100.000
recruitment.YCS116   0.742500.742500.000
recruitment.YCS117   0.570500.570500.000
recruitment.YCS118   0.451600.451600.000
recruitment.YCS119   1.000001.000000.000
recruitment.YCS120   1.000001.000000.000
size_at_age.Linf_male47.1200047.120000.000
size_at_age.cv2_male 0.062560.062560.000
selectivity[potSsel].male1    3.039003.039000.000
selectivity[potSsel].male2    0.190700.190600.052
selectivity[potSsel].female1  3.316003.316000.000
selectivity[potSsel].female2  1.045001.045000.000
Casal2 parameter estimates
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
process[Recruitment].b0              13871.0000013879.1000013871.0000013871.1000013879.1000013871.00000
age_length[asMm0].linf            46.6612046.6574046.6613046.6613046.6573046.66130
age_length[asMm0].cv_last         0.020270.022600.020270.020280.022600.02027
selectivity[potSurveySel_male].a50             3.043133.039273.038523.025473.048233.02520
selectivity[potSurveySel_male].ato95           0.213190.191340.190390.125870.234970.12457
selectivity[potSurveySel_female].a50             3.324033.324883.324033.324023.324893.32403
selectivity[potSurveySel_female].ato95           1.078861.072241.078871.078871.072251.07887
process[Recruitment].ycs_values{1980}0.212020.211930.212020.211990.211900.21202
process[Recruitment].ycs_values{1981}0.208850.208780.208850.208760.208770.20885
process[Recruitment].ycs_values{1982}0.209410.209390.209410.209540.209420.20941
process[Recruitment].ycs_values{1983}0.217600.217620.217600.217600.217630.21760
process[Recruitment].ycs_values{1984}0.240420.240460.240420.240420.240460.24042
process[Recruitment].ycs_values{1985}0.292900.292900.292910.292910.292900.29291
process[Recruitment].ycs_values{1986}0.400430.400290.400430.400430.400290.40043
process[Recruitment].ycs_values{1987}0.600380.600050.600390.600380.600040.60039
process[Recruitment].ycs_values{1988}0.926670.926230.926610.926600.926310.92661
process[Recruitment].ycs_values{1989}1.475261.475321.475211.475221.475131.47521
process[Recruitment].ycs_values{1990}1.596141.596971.596251.596251.597181.59625
process[Recruitment].ycs_values{1991}1.004721.005311.004681.004701.005251.00468
process[Recruitment].ycs_values{1992}0.835520.835940.835560.835560.835940.83556
process[Recruitment].ycs_values{1993}0.852480.852510.852510.852500.852510.85251
process[Recruitment].ycs_values{1994}1.110681.109491.110521.110531.109451.11052
process[Recruitment].ycs_values{1995}1.664561.663181.664581.664561.663321.66458
process[Recruitment].ycs_values{1996}1.053271.054041.053341.053361.053951.05334
process[Recruitment].ycs_values{1997}0.823920.823980.823950.823940.823990.82395
process[Recruitment].ycs_values{1998}0.843940.843900.843900.843890.843910.84390
process[Recruitment].ycs_values{1999}0.621690.621890.621680.621680.621890.62168
process[Recruitment].ycs_values{2000}0.741840.742360.741840.741840.742360.74184
process[Recruitment].ycs_values{2001}0.673330.673560.673330.673330.673560.67333
process[Recruitment].ycs_values{2002}0.443960.443960.443960.443960.443960.44396
process[Recruitment].ycs_values{2003}1.210101.207401.210101.210091.207401.21010
process[Recruitment].ycs_values{2004}1.220061.215571.220051.220051.215571.22005
process[Recruitment].ycs_values{2005}0.489320.487250.489320.489320.487250.48932
process[Recruitment].ycs_values{2006}1.003641.000001.003631.003631.000001.00363
process[Recruitment].ycs_values{2007}0.711870.722990.711870.711870.723000.71187
process[Recruitment].ycs_values{2008}0.831320.830480.831320.831320.830480.83132
process[Recruitment].ycs_values{2009}0.421530.421280.421530.421530.421280.42153
process[Recruitment].ycs_values{2010}1.261311.260961.261311.261311.260961.26131
process[Recruitment].ycs_values{2011}0.665270.666030.665270.665270.666030.66527
process[Recruitment].ycs_values{2012}0.709380.709520.709380.709380.709520.70938
process[Recruitment].ycs_values{2013}0.551840.551650.551840.551840.551650.55184
process[Recruitment].ycs_values{2014}0.429670.429380.429670.429670.429380.42967
Casal2 parameter estimates: Percent Difference from betadiff_casal_flags_on
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
process[Recruitment].b0              0.000-0.0580.000-0.001-0.0580.000
age_length[asMm0].linf            0.0000.008-0.000-0.0000.008-0.000
age_length[asMm0].cv_last         0.000-11.507-0.001-0.079-11.481-0.001
selectivity[potSurveySel_male].a50             0.0000.1270.1510.580-0.1680.589
selectivity[potSurveySel_male].ato95           0.00010.24910.69840.957-10.21441.569
selectivity[potSurveySel_female].a50             0.000-0.0260.0000.000-0.0260.000
selectivity[potSurveySel_female].ato95           0.0000.614-0.001-0.0010.613-0.001
process[Recruitment].ycs_values{1980}0.0000.0400.0000.0120.0540.000
process[Recruitment].ycs_values{1981}0.0000.0330.0000.0440.0370.000
process[Recruitment].ycs_values{1982}0.0000.0120.000-0.058-0.0020.000
process[Recruitment].ycs_values{1983}0.000-0.0100.000-0.000-0.0120.000
process[Recruitment].ycs_values{1984}0.000-0.016-0.000-0.000-0.016-0.000
process[Recruitment].ycs_values{1985}0.0000.002-0.001-0.0010.002-0.001
process[Recruitment].ycs_values{1986}0.0000.034-0.0000.0000.034-0.000
process[Recruitment].ycs_values{1987}0.0000.054-0.002-0.0010.055-0.002
process[Recruitment].ycs_values{1988}0.0000.0480.0060.0070.0390.006
process[Recruitment].ycs_values{1989}0.000-0.0040.0030.0030.0090.003
process[Recruitment].ycs_values{1990}0.000-0.052-0.007-0.007-0.065-0.007
process[Recruitment].ycs_values{1991}0.000-0.0590.0040.002-0.0530.004
process[Recruitment].ycs_values{1992}0.000-0.051-0.005-0.005-0.051-0.005
process[Recruitment].ycs_values{1993}0.000-0.004-0.003-0.003-0.004-0.003
process[Recruitment].ycs_values{1994}0.0000.1070.0140.0140.1110.014
process[Recruitment].ycs_values{1995}0.0000.083-0.0010.0000.074-0.001
process[Recruitment].ycs_values{1996}0.000-0.073-0.007-0.009-0.065-0.007
process[Recruitment].ycs_values{1997}0.000-0.008-0.004-0.003-0.009-0.004
process[Recruitment].ycs_values{1998}0.0000.0040.0040.0050.0030.004
process[Recruitment].ycs_values{1999}0.000-0.0320.0010.001-0.0330.001
process[Recruitment].ycs_values{2000}0.000-0.0710.000-0.000-0.0710.000
process[Recruitment].ycs_values{2001}0.000-0.0340.0010.000-0.0340.001
process[Recruitment].ycs_values{2002}0.0000.0010.0000.0000.0010.000
process[Recruitment].ycs_values{2003}0.0000.2230.0000.0010.2230.000
process[Recruitment].ycs_values{2004}0.0000.3680.0010.0010.3680.001
process[Recruitment].ycs_values{2005}0.0000.4210.0000.0000.4210.000
process[Recruitment].ycs_values{2006}0.0000.3630.0010.0010.3630.001
process[Recruitment].ycs_values{2007}0.000-1.5630.0000.000-1.5630.000
process[Recruitment].ycs_values{2008}0.0000.1010.0000.0000.1010.000
process[Recruitment].ycs_values{2009}0.0000.0590.0000.0000.0590.000
process[Recruitment].ycs_values{2010}0.0000.0280.0000.0000.0280.000
process[Recruitment].ycs_values{2011}0.000-0.1140.000-0.000-0.1140.000
process[Recruitment].ycs_values{2012}0.000-0.0190.000-0.000-0.0190.000
process[Recruitment].ycs_values{2013}0.0000.0350.0000.0000.0340.000
process[Recruitment].ycs_values{2014}0.0000.0670.000-0.0000.0670.000

CASAL objective function component values
ComponentBase_ModelSensitivity_1
potCPUE-51.180-51.180
potSurveyAbd-3.254-3.254
logbookLF54.66054.660
potCommAF54.41054.410
potSurveyAF118.200118.200
recCatchLF68.04068.040
prior_on_initialization.B09.5599.559
prior_on_recruitment.YCS9.6589.658
prior_on_size_at_age.Linf_male0.0000.000
prior_on_size_at_age.cv2_male0.0000.000
prior_on_selectivity[potSsel].male0.0000.000
prior_on_selectivity[potSsel].female0.0000.000
prior_on_q_qPotCPUE-8.331-8.331
prior_on_q_potSurveyq-11.380-11.380
clp1F0.0000.000
clp2F0.0000.000
clp3F0.0000.000
YCS0.5660.566
Total241.000241.000
Casal2 objective function component values
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_toladolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
observation->potSurvey-2010                 -0.789-0.793-0.789-0.789-0.793-0.789
observation->potSurvey-2014                 -0.931-0.933-0.931-0.931-0.933-0.931
observation->potSurvey-2018                 -1.120-1.121-1.120-1.120-1.121-1.120
observation->potCPUE-1990                 1.1261.1291.1261.1251.1291.126
observation->potCPUE-1991                 2.2272.2332.2272.2272.2332.227
observation->potCPUE-1992                 -0.686-0.682-0.686-0.686-0.682-0.686
observation->potCPUE-1993                 -1.938-1.937-1.938-1.938-1.937-1.938
observation->potCPUE-1994                 -2.140-2.140-2.140-2.140-2.140-2.140
observation->potCPUE-1995                 -2.056-2.056-2.056-2.056-2.056-2.056
observation->potCPUE-1996                 -2.194-2.194-2.194-2.194-2.194-2.194
observation->potCPUE-1997                 -2.026-2.025-2.026-2.026-2.025-2.026
observation->potCPUE-1998                 -2.164-2.165-2.164-2.164-2.165-2.164
observation->potCPUE-1999                 -1.736-1.734-1.736-1.736-1.734-1.736
observation->potCPUE-2000                 -2.194-2.194-2.194-2.194-2.194-2.194
observation->potCPUE-2001                 -2.173-2.173-2.173-2.173-2.173-2.173
observation->potCPUE-2002                 -2.084-2.084-2.084-2.084-2.084-2.084
observation->potCPUE-2003                 -2.187-2.187-2.187-2.187-2.187-2.187
observation->potCPUE-2004                 -2.185-2.185-2.185-2.185-2.185-2.185
observation->potCPUE-2005                 -1.894-1.895-1.894-1.894-1.895-1.894
observation->potCPUE-2006                 -1.596-1.601-1.596-1.596-1.601-1.596
observation->potCPUE-2007                 -1.999-2.003-1.999-1.999-2.003-1.999
observation->potCPUE-2008                 -1.932-1.936-1.932-1.932-1.936-1.932
observation->potCPUE-2009                 -1.631-1.633-1.631-1.631-1.633-1.631
observation->potCPUE-2010                 -2.186-2.186-2.186-2.186-2.186-2.186
observation->potCPUE-2011                 -1.876-1.866-1.876-1.876-1.866-1.876
observation->potCPUE-2012                 -1.908-1.898-1.908-1.908-1.898-1.908
observation->potCPUE-2013                 -1.959-1.955-1.959-1.959-1.955-1.959
observation->potCPUE-2014                 -1.313-1.311-1.313-1.313-1.311-1.313
observation->potCPUE-2015                 -1.908-1.909-1.908-1.908-1.909-1.908
observation->potCPUE-2016                 -2.082-2.083-2.082-2.082-2.083-2.082
observation->potCPUE-2017                 -2.037-2.039-2.037-2.037-2.039-2.037
observation->potCPUE-2018                 -1.654-1.650-1.654-1.654-1.650-1.654
observation->AFpotSurvey-2010                 35.03635.08535.03635.03635.08535.036
observation->AFpotSurvey-2014                 38.71838.72038.71838.71838.72038.718
observation->AFpotSurvey-2018                 38.82138.79138.82138.82138.79138.821
observation->AFpotFishing-2018                 45.61145.65345.61145.61145.65345.611
observation->AFpotFishing-2019                 7.8547.8547.8547.8547.8547.854
observation->lgobookLF                 51.53851.54151.53851.53751.54151.538
observation->rec.catch.LF              71.22272.11371.22271.22272.11371.222
prior->B0->process[Recruitment].b0              9.5389.5389.5389.5389.5389.538
prior->male_Linf->age_length[asMm0].linf            0.0000.0000.0000.0000.0000.000
prior->male_cv2->age_length[asMm0].cv_last         0.0000.0000.0000.0000.0000.000
prior->potSurvey_mA50->selectivity[potSurveySel_male].a50             0.0000.0000.0000.0000.0000.000
prior->potSurvey_mAto95->selectivity[potSurveySel_male].ato95           0.0000.0000.0000.0000.0000.000
prior->potSurvey_fA50->selectivity[potSurveySel_female].a50             0.0000.0000.0000.0000.0000.000
prior->potSurvey_fAto95->selectivity[potSurveySel_female].ato95           0.0000.0000.0000.0000.0000.000
prior->YCS->process[Recruitment].ycs_values{1980}1.6241.6251.6241.6241.6261.624
prior->YCS->process[Recruitment].ycs_values{1981}1.6781.6791.6781.6791.6791.678
prior->YCS->process[Recruitment].ycs_values{1982}1.6681.6681.6681.6661.6681.668
prior->YCS->process[Recruitment].ycs_values{1983}1.5331.5331.5331.5331.5331.533
prior->YCS->process[Recruitment].ycs_values{1984}1.2041.2041.2041.2041.2041.204
prior->YCS->process[Recruitment].ycs_values{1985}0.6480.6480.6480.6480.6480.648
prior->YCS->process[Recruitment].ycs_values{1986}0.0280.0280.0280.0280.0280.028
prior->YCS->process[Recruitment].ycs_values{1987}-0.304-0.303-0.304-0.304-0.303-0.304
prior->YCS->process[Recruitment].ycs_values{1988}-0.066-0.067-0.066-0.066-0.067-0.066
prior->YCS->process[Recruitment].ycs_values{1989}0.8680.8680.8670.8670.8670.867
prior->YCS->process[Recruitment].ycs_values{1990}1.0951.0971.0961.0961.0971.096
prior->YCS->process[Recruitment].ycs_values{1991}0.0460.0460.0450.0460.0460.045
prior->YCS->process[Recruitment].ycs_values{1992}-0.179-0.178-0.179-0.179-0.178-0.179
prior->YCS->process[Recruitment].ycs_values{1993}-0.160-0.160-0.160-0.160-0.160-0.160
prior->YCS->process[Recruitment].ycs_values{1994}0.2140.2120.2140.2140.2120.214
prior->YCS->process[Recruitment].ycs_values{1995}1.2251.2221.2251.2251.2231.225
prior->YCS->process[Recruitment].ycs_values{1996}0.1210.1220.1210.1210.1220.121
prior->YCS->process[Recruitment].ycs_values{1997}-0.191-0.191-0.191-0.191-0.191-0.191
prior->YCS->process[Recruitment].ycs_values{1998}-0.169-0.169-0.169-0.169-0.169-0.169
prior->YCS->process[Recruitment].ycs_values{1999}-0.307-0.307-0.307-0.307-0.307-0.307
prior->YCS->process[Recruitment].ycs_values{2000}-0.264-0.264-0.264-0.264-0.264-0.264
prior->YCS->process[Recruitment].ycs_values{2001}-0.300-0.300-0.300-0.300-0.300-0.300
prior->YCS->process[Recruitment].ycs_values{2002}-0.107-0.107-0.107-0.107-0.107-0.107
prior->YCS->process[Recruitment].ycs_values{2003}0.3840.3790.3840.3840.3790.384
prior->YCS->process[Recruitment].ycs_values{2004}0.4010.3930.4010.4010.3930.401
prior->YCS->process[Recruitment].ycs_values{2005}-0.203-0.199-0.203-0.203-0.199-0.203
prior->YCS->process[Recruitment].ycs_values{2006}0.0440.0380.0440.0440.0380.044
prior->YCS->process[Recruitment].ycs_values{2007}-0.284-0.277-0.284-0.284-0.277-0.284
prior->YCS->process[Recruitment].ycs_values{2008}-0.183-0.184-0.183-0.183-0.184-0.183
prior->YCS->process[Recruitment].ycs_values{2009}-0.044-0.043-0.044-0.044-0.043-0.044
prior->YCS->process[Recruitment].ycs_values{2010}0.4740.4740.4740.4740.4740.474
prior->YCS->process[Recruitment].ycs_values{2011}-0.303-0.303-0.303-0.303-0.303-0.303
prior->YCS->process[Recruitment].ycs_values{2012}-0.285-0.285-0.285-0.285-0.285-0.285
prior->YCS->process[Recruitment].ycs_values{2013}-0.279-0.278-0.279-0.279-0.278-0.279
prior->YCS->process[Recruitment].ycs_values{2014}-0.068-0.067-0.068-0.068-0.067-0.068
additional_prior->PrSurveyQ                 -11.153-11.154-11.153-11.153-11.154-11.153
additional_prior->PrCPUEQ                 -8.207-8.208-8.207-8.207-8.208-8.207
additional_prior->YCS_average_1                 0.4950.4960.4950.4950.4960.495
total_negloglike                 237.802238.776237.802237.802238.776237.802
total_score                 237.802238.776237.802237.802238.776237.802

## [1] "CASAL base model convergence information"
## [1] "Minimiser convergence threshold: 2e-005"                                                                 
## [2] "fmm: have converged: t = 1.53773e-005 f = 240.986"                                                       
## [3] "Successful convergence in optimise"                                                                      
## [4] "Minimiser achieved convergence after 91 quasi-Newton iterations using 101 objective function evaluations"
## [1] ""
## [1] "CASAL sensitivity 1 model convergence information"
## [1] "Minimiser convergence threshold: 1e-007"                                                                                 
## [2] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [3] "Successful convergence in optimise"                                                                                      
## [4] "Minimiser achieved convergence after 99 quasi-Newton iterations using 119 objective function evaluations"                
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model convergence information"
## [1] "fmm: have converged: t = 1.80811e-05 f = 237.802"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model convergence information"
## [1] "fmm: have converged: t = 1.05483e-05 f = 238.776"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model convergence information"
## [1] "fmm: have converged: t = 2.57984e-09 f = 237.802"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model convergence information"
## character(0)
## [1] "CASAL model warnings"
## NULL
## [1] ""
## [1] "CASAL sensitivity 1 model warnings"
## NULL
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [8] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [9] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [8] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [9] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [8] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [9] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model warnings"
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"                                                                              
## [2] "Cannot report the length-based selectivity values. This report (__lineFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"     
## [3] "Cannot report the length-based selectivity values. This report (__lineFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"
## [4] "Cannot report the length-based selectivity values. This report (__potFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [5] "Cannot report the length-based selectivity values. This report (__potFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [6] "Cannot report the length-based selectivity values. This report (__recFretained__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information" 
## [7] "Cannot report the length-based selectivity values. This report (__recFSel__) is being ignored. This can be done using the Casal2 R package. See the User Manual for more information"      
## [8] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"                                                                                                      
## [9] "estimated parameter 'age_length[asMm0].cv_last' was within 0.001 of lower bound 0.02"

1.5 Matching of outputs

Time series comparisons with CASAL base model results

## [1] "Catch time series base model comparison for run betadiff_casal_flags_on"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_off"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_on_low_tol"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_off"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on_low_tol"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""

Derived quantities

SB0, SBcurrent, MSY, F_MSY, others…

1.6 Plots

Comparison plots

## [1] "CASAL base parameter correlation range (excluding 1.0): -0.6876 0.9999"

## [1] "CASAL sensitivity 1 parameter correlation range (excluding 1.0): -0.6894 0.9999"

## [1] "Casal2 betadiff_casal_flags_on parameter correlation range (excluding 1.0): -0.704682058267449 0.999819516983665"

## [1] "Casal2 betadiff_casal_flags_off parameter correlation range (excluding 1.0): -0.697745743037247 0.999863100354264"

## [1] "Casal2 betadiff_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.722472902754433 0.999845201096008"

## [1] "Casal2 adolc_casal_flags_on parameter correlation range (excluding 1.0): -0.704683609734695 0.99981954195823"

## [1] "Casal2 adolc_casal_flags_off parameter correlation range (excluding 1.0): -0.697745743037247 0.999863100354264"

## [1] "Casal2 adolc_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.719263523071944 0.999846316745602"